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We study the localization properties in coupled double quantum wells with an in-plane magnetic 
field. The localization length is directly calculated using a transfer matrix technique and finite size 
scaling analysis. We show that the system maps into a 2D XY model and undergoes a disorder driven 
Kosterlitz-Thouless type metal-insulator transition depending on the coupling strength between the 
two-dimensional layers and the magnitude of the in-plane magnetic field. For a system with fixed 
disorder, the metallic regime appears to be a window in the magnetic field - coupling strength 
plane. Experimental implications of the transition will be discussed. 
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There have been long lasting interests to understand 
localization problem in two-dimensional(2D) systems [Q. 
According to the scaling theory of localization 0, all 
states in a two-dimensional system are localized if only 
scalar random potential is present. However, in the pres- 
ence of a strong perpendicular magnetic field, extended 
states appear in the center of disorder broadened Landau 
bands and give arise to the quantum Hall effect || . There 
has been report on metal-insulator transition (MIT) in 
the presence of spin-orbit coupling in two-dimensional 
systems B . It is important to understand the criteria for 
MIT in disordered 2D systems and corresponding univer- 
sality classes for these phase transitions. 

From symmetry consideration, 2D localization prob- 
lem can be classified into following universality classes: 
(1) orthogonal ensemble, where only scalar random po- 
tential is present; (2) symplectic ensemble, such as the 
case with strong spin-orbit scattering; (3) unitary ensem- 
ble, such as the quantum Hall system, where the strong 
perpendicular magnetic field breaks the time-reversal 
symmetry. One interesting system belongs to the uni- 
tary class is the system with random mag netic field §. 
This system has attracted much attention recently be- 
cause it is relevant to the studies of fractional quantum 
Hall systems at even denominator filling factors and high- 
T c models with gauge field fluctuations. In this prob- 
lem, the time-reversal symmetry is locally broken by the 
random magnetic field, whether this symmetry-breaking 
is sufficient to induce metal-insulator transition remains 
unclear despite extensive recent studies |5J. However, 
Zhang and Arovas [J7J developed an effective field theory 
for the random field problem, where they found the MIT 
in 2D random field is a Kosterlitz-Thouless type transi- 
tion§. 

In this letter, we study a different system where the 
time-reversal symmetry is also broken: coupled double 
quantum wells with in-plane magnetic field. This sys- 



tem is experimentally more accessible comparing with the 
random field system and theoretical predictions can be 
tested. We will demonstrate that the system maps into a 
2D random field problem and a Kosterlitz-Thouless (KT) 
type metal-insulator transition arises depending on disor- 
der, the coupling strength between the two-dimensional 
layers and the magnitude of the in-plane magnetic field. 
With fixed disorder, the metallic regime will be a window 
in the magnetic field - coupling strength plane. Experi- 
mental implications of the transition will be discussed. 




FIG. 1. Sketch of double quantum wells with in-plane mag- 
netic field. 

Consider an electron travels through the system of cou- 
pled double quantum wells with in-plane magnetic field. 
The magnetic field is only relevant if the electron path 
contains closed loops through inter-layer hopping (Loops 
3 and 4 in Fig.l.). However, the angles between these 
loops and the in-plane magnetic field are random, so the 
relevant magnetic field components along the direction 
perpendicular to these paths are random and its sign can 
be either positive or negative. In another words, this 
system topologically maps into a 2D random field prob- 
lem. Those paths staying within each layer (loops 1 and 
2 in Fig.l) or not containing closed loops correspond to 
regions with zero field in the random field problem. 
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According to the effective field theory put forward by 
Zhang and Arovas JtJ for the 2D random field problem, 
a long range logarithmic interaction arises between the 
topological density of the nonlinear a model, or more 
physically the edge excitations of the fluctuating mag- 
netic domains in the random magnetic field case; and in 
the coupled double quantum well system, they should be 
the hopping paths with different orientations regarding to 
the in-plane magnetic field. They demonstrated that the 
2D random field problem at zero temperature could map 
into a 2D XY model at finite temperature [Q in which 
the longitudinal conductance corresponds to the inverse 
of the temperature. As opposed to the 2D XY model 
there should be a disorder driven (where strong disorder 
corresponds to small conductance, hence high tempera- 
ture) Kosterlitz-Thouless transition in 2D random field 
problem, or in our case, the system of coupled double 
quantum wells with in-plane magnetic field at zero tem- 
perature. Direct consequence of such transition is a line 
of fixed points with continuously varying conductance in 
the metallic regime . Using a transfer matrix technique 
and finite-size scaling analysis, we directly calculate the 
localization length in the system to demonstrate the ex- 
istence of the line of fixed points and critical behavior 
consistent with KT transition. 

In the following, we briefly outline our model and the 
technique we used to calculate the finite-size localization 
lcngth(AM) and to extrapolate the thermodynamic local- 
ization length(£). We model our bi-layer system as two 
square lattice with interlayer hopping t p . The square lat- 
tice is a very long strip geometry with a finite width (M) 
square lattice with nearest neighbor hopping. Periodic 
boundary condition in the width direction is assumed. 
The disorder potential is modeled by the on-site white- 
noise potential Vi rn i (i denotes the column index, m de- 
notes the chain index, I is the layer index) ranging from 
— W/2 to W/2. The effect of the in-plane magnetic field 
appears in the complex phase of the inter-layer hopping 
term. The strength of the magnetic field is characterized 
by the flux per plaquette (</>) in unit of magnetic flux 
quanta (<f) = hc/e). The Hamiltonian of this system can 
be written as: 

M 

U = X) X] v imi\iml >< iml\ (1) 

2 = 1,2 i m=l 

-t V ^2 [\iml >< jnl\ + \jnl >< iml\] 

1—1,2 <im;jn> 

- \^ [tp(i)\iml X im2\ + t P (i)\im2 >< iml\] , 

im 

where < im; jn > indicates nearest neighbors on the lat- 
tice. The amplitude of the in-layer hopping term(i) is 
chosen as the unit of the energy. For a specific energy 
E, a transfer matrix Ti can be easily set up mapping the 
wavefunction amplitudes at column i — 1 and i to those 
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FIG. 2. (a) Xm/M as function of disorder strength W 
without magnetic field for M — 4 (circles), M — 8 (stars), 
M = 12 (diamonds), M = 16 (triangles) and M — 24 
(squares). Inset: Thermodynamic localization length £ as 
function of disorder W. (b) Xm/M as function of disorder 
strength W with magnetic field (f> — 1/5 for M — 4 (circles), 
M = 8 (stars), M = 12 (diamonds), M = 16 (triangles) and 
M — 24 (squares). Inset: Thermodynamic localization length 
£ as function of disorder W . (c) Fit the thermodynamic local- 
ization length £ with the KT transition and power-law expres- 
sions. Inset: scaling function using the data from Fig. 2(b). 

at column Using a standard iteration algorithm [|l0| , 
we can calculate the Lyapunov exponents for the trans- 
fer matrix X^. The localization length Xm(E) for energy 
E at finite width M is then given by the inverse of the 
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smallest Lyapunov exponent. In our numerical calcula- 
tion, we choose the sample length to be over 10 5 so that 
the self-averaging effect automatically takes care of the 
ensemble statistical fluctuations. We use the standard 
one-parameter finite-size scaling analysis |i"l[ to obtain 
the thermodynamic localization length £. The scaling 
calculations are carried out with varying disorder W (or 
in-plane magnetic field <f> 1 or inter-layer hopping ampli- 
tude t p ) at fixed energy. 

We now discuss our numerical results in various situa- 
tions and their implications. Presented in Fig. 2(a) is the 
ratio of finite localization length Am with system width 
M versus disorder W without in-plane magnetic field. 
Energy is set at the band center with E = 0. Different 
curves are for different system widths(M=4-24). Within 
numerical fluctuations, the curve for smaller M is above 
that of the larger M for the entire range of W. This 
indicates that there is no metal-insulator transition for 
finite disorder W, consistent with one's expectation for 
the double-layer system with zero magnetic field. In the 
inset, we show the thermodynamic localization length £ 
as function of disorder strength W. Although £ rises fast 
for small W, numerical fitting still shows that the critical 
value W c ~ 0. 

Fig. 2(b) shows the main results of our paper. Here we 
plot the same curves as in Fig. 2(a) but with finite in- 
plane field 4> = 1/5. The striking difference with Fig. 2(a) 
is that all curves merge together for W < W c ~ 5.5. As 
is well known from the finite-size scaling studies of phase 
transition, all curves for different sizes should cross at 
single point (critical point) for a conventional continuous 



Fig. 3(b), we see all curves merge together in the range 
t p ~ 0.1 — 3.0, indicating that we have continuous critical 
points in this range. 



transition |12|. Two examples for such transition in the 
localization problem are the three-dimensional Anderson 
model p] ] and two-dimensional case with spin-orbital 
interaction However, Fig. 2(b) is quite different from 
the conventional transition in the sense that there is no 
single crossing point, but all curves merge together for 
W < W c . This shows that all the points for W < W c 
are critical points |]l2| , namely, we have a line of critical 
points which is consistent with our argument that the 
system undergoes a disorder driven Kosterlitz-Thouless 
transition. In Fig. 2(c), we fit the data for £ with £ cx 
exp(a/y/W — W c ) (full circles and solid line), typical for 
the KT transition, and £ cx (W — W c ) u (open circles and 
dotted line) for W > W c . For the KT fitting, W c = 5.7 
and a = 13.8 ± 0.1 and for the power-law fitting, W c = 
6.2 and v = 3.4 ± 0.3. Results clearly show that KT 
dependence has a better fitting and supports the notion 
that the transition we found in Fig. 2(b) belongs to KT 
transition. 

In Fig. 3 we show further numerical support for KT 
transition by ploting Xm/M versus interlayer hopping 
t p for fixed disorder. Fig. 3(a) is for zero magnetic field 
and we see that curve with smaller M has larger value 
of Xm/M for all t p , consistent with the picture that all 
states are localized. For finite field with <f> — 1/5 in 
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FIG. 3. (a) Xm/M as function of inter-layer hopping am- 
plitude t p without magnetic field for M = 8 (circles) , M = 12 
(stars), M — 16 (diamonds), M = 24 (triangles), (b) X M /M 
as function of inter-layer hopping amplitude t v with magnetic 
field <j> = 1/5 for M = 8 (circles), M = 12 (stars), M = 16 
(diamonds), M = 24 (triangles). Inset: Phase diagram on the 
4> — t p plane. The filled circles are for W=6, the opened ones 
are for W=5. The states within each closed loop are extended 
states. 

In a realistic experimental set-up, the sample disor- 
der is fixed and the in-plane magnetic field can be easily 
changed. In the following we study the phase diagram 
for fixed disorder by varying magnetic field and inter- 
layer hopping which depends on the distance between 
the two layers. For a fixed moderate coupling strength 
between two layers, in-plane magnetic has no effect on 
those electron paths within each layer, these states tend 
to stay localized. However, for those paths going through 
inter-layer hopping, the time-reversal symmetry is bro- 
ken. So as the magnitude of magnetic field increases from 
zero, the system undergoes a metal-insulator transition 
at a critical magnetic field. In the presence of extremely 
strong magnetic field, the magnetic length is compara- 
ble to the inter-layer separation, the system will become 



3 



localized again with all electrons quenched within the 
magnetic length. The story is similar in a fixed moder- 
ate magnetic field. When the inter-layer coupling is very 
small, the number of paths going through inter-layer hop- 
ping is very limited. So the system behaves as two sep- 
arate layers and all states are localized. When the cou- 
pling strength is very large, the energy gap between the 
symmetric and antisymmetric states is very large. Thus, 
the hopping between the symmetric and antisymmetric 
states is very weak, and again all states are localized. Ex- 
tended states only emerge in the intermediate coupling 
regime. We have carried out numerical calculations for 
many points on the <j) — t p plane. Presented in the insert 
of Fig. 3(b) is the resulting phase boundary separating 
the metallic and insulator regime. The phase diagram is 
determined by the points when those lines of Xm/M for 
different M merge together. Indeed the metallic regime 
is a window on the <f> — t p plane and it shrinks when the 
system gets more disordered. 



E = 0.0, W=7.0, t=2.0 
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FIG. 4. Average conductance (over 100 samples) G as 
function of sample size. Different symbol represents different 
in-plane magnetic field as shown in the inset (field dependent 
conductance for M = 32). 

The above phase diagram can be tested by measuring 
the temperature dependence of the conductance in cou- 
pled double quantum wells with in-plane magnetic field. 
We mimic this effect by calculating the sample size de- 
pendence of such conductance, since finite temperature 
provides an effective cut-off length (the inelastic scatter- 
ing length) in a real system. The conductance is calcu- 
lated using the Landauer-Buttiker formula and recursive 
Green's function technique. As presented in Fig.^, one 
can see that the conductance in the metallic regime in- 
creases with the sample size(mimic decreasing tempera- 
ture), while that for insulator regime decreases. The con- 
ductance should reach maximum in the extended regime 
while varying the in-plane field at fixed sample size (fixed 
temperature) as presented in the inset. To test the ef- 
fect of the coupling strength, one could do similar ex- 
periment on different samples with various barrier height 



and thickness in the double quantum wells. 

In summary, we have demonstrated that the system of 
double quantum wells with in-plane magnetic field maps 
into the 2D random field problem and hence a 2D XY 
model, which undergoes a disorder driven KT type metal- 
insulator transition depending the inter-layer coupling 
and the magnetic field. In fact, the metallic regime is 
a window on the field-coupling plane which shrinks with 
increasing disorder. We propose temperature dependence 
transport measurement to test this phase diagram. 

We thank Song He and Michael Ma for many helpful 
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